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Symmetry of the crystal lattice can be a determining factor for the structure of Cooper pairs 
in unconventional superconductors. In this study we extend the discussion of superconductivity in 
non-centrosymmetric materials to the case when inversion symmetry is missing locally, but is present 
on a global level. Concretely, we investigate the staggered non-centrosymmetricity within a regular 
sublattice structure, in some analogy to the discussion of superconductivity in antiferromagnetic 
systems. Three crystal structures are analyzed in detail as illustrative examples for the extended 
classification of Cooper-pairing channels. One of the cases may be relevant for the class of iron- 
pnictide superconductors. 
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I. INTRODUCTION 

Shortly after the seminal paper by Bardeen, Cooper, 
and Schrieffer describing superconductivity through pair- 
ing of electrons of equal energy and opposite spin and 
momentum,— Anderson realized that the existence of 
such degenerate electron pairs would be guaranteed quite 
generally by time-reversal symmetry^ Indeed, remov- 
ing time-reversal symmetry by an external magnetic 
field, magnetic impurities, or ferromagnetic order sub- 
stantially weakens or even suppresses superconductivity 
in the spin-singlet channel. Later, Baltensperger and 
Strassler demonstrated that spin-singlet superconductiv- 
ity and antifcrromagnctism can coexist for an appropri- 
ate pair structured In such systems, staggered moments 
break time-reversal symmetry only on sublattices. How- 
ever, time-reversal operation may be undone globally by 
exchanging the two sublattices. In this case, the spin- 
singlet pair wavefunction has dominant amplitudes for 
the two electrons being on different sublattices. 

For Cooper pairing of electrons in the spin-triplet 
configuration, Anderson showed several years later that 
an additional discrete symmetry is needed, namely in- 
version symmetry^ The discovery of superconductivity 
in crystals lacking an inversion center and yet show- 
ing features usually attributed to spin-triplet pairing, 
therefore, attracted much attention in recent years. 
Non-centrosymmetricity affects the electronic spectrum 
through symmetry specific antisymmetric spin-orbit cou- 
pling. Spin-triplet superconductivity is not simply sup- 
pressed in favor of spin-singlet pairing, but actually elec- 
trons pair with a mixed-parity structure combining a 
spin-singlet and -triplet component. 

Since in the context of time-reversal-symmetry break- 
ing and superconductivity the effect of both, ferromag- 
netic and antiferromagnetic order on the Cooper-pair for- 
mation have been studied, it is natural to extend the 
recent discussion of globally non-centrosymmetric super- 
conductivity to its staggered form. Recently, Yanase has 



analyzed the case of locally broken inversion symmetry 
due to stacking faults, where the global inversion sym- 
metry is retained because of the random distribution of 
these faults^ In the present study, we will generalize the 
discussion of broken inversion symmetry from a ferro- 
type to the antiferro-type. 

After introducing first a general formulation for the 
antiferro-type of "non-centrosymmetric" lattices, we will 
discuss two examples in detail to illustrate the influence 
on superconductivity and then apply the results to a 
crystal structure as found in the iron pnictide supercon- 
ductors. The underlying crystal symmetry for all three 
systems is tetragonal and can be characterized by a spe- 
cific sublattice structure of two distinct types of sites or 
bonds yielding a doubling of the ordinary unit cell. Each 
of them has a different subgroup of D^h leaving the sub- 
lattice structure invariant. The examples then differ in 
that the first and third example have a sublattice lacking 
inversion symmetry, while the second example lacks in- 
version symmetry only on the bonds connecting the two 
sublattices. 



II. SINGLE PARTICLE HAMILTONIAN 

Before looking at these specific examples, we introduce 
a general formalism for superconductivity in a lattice 
with a non-centrosymmetric sublattice structure. While 
such a crystal has centers of inversion, the lattice struc- 
ture includes local violations of inversion symmetry (see 
Fig. [T|), which yield a staggered form of antisymmetric 
spin-orbit coupling. This can be incorporated into the 
kinetic energy by defining a folded Brillouin zone with 
two bands characterized by the wave vector Q (2Q is a 
reciprocal lattice vector). Thus, we define the operators, 



Cks a = 1 
Ck+Qs a = 2, 



(1) 



where we use a — 1, 2 as band indices. 
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FIG. 1: The two example crystal structures analyzed in de- 
tail in Sees. Hil l and IIVI (a) Inversion-symmetry lacking layers 
that are stacked along the z direction in a staggered way. The 
symmetry center lies between the layers and the crystal has 
a symmorphic structure, (b) Top view of the crystal struc- 
ture with 06 octahedra rotated around the c axis leading to 
a doubling of the unit cell. While both sublattices still retain 
inversion symmetry, the bonds do not as the rotation shifts 
the O ions off the bonds. This crystal structure with its sym- 
metry center on one of the sublattices is non-symmorphic. 



A. Diagonal single-particle Hamiltonian 



in the Hamiltonian of the form 

H = ^A M (c t lks c 2ks + h.c.) (5) 

k,s 

or again in terms of Pauli matrices 

E w = A/iT^,, (6) 

off-diagonal in the band index. Note that time-reversal 
symmetry leads to the condition that 

^koa' — ^ — kaa' (7) 

and Skaa' is spin independent. 

The four Pauli matrices (r° , f) for the band part are 
easily interpreted in the sublattice and two-band notion. 
Matrix elements independent of sublattice arc diagonal 
represented by r° for intra-sublattice and by r 3 for inter- 
sublattice processes. Analogously, inter-band hybridiza- 
tion is incorporated in r 1 and r 2 for intra- and inter- 
sublatticc couplings, respectively. This is summarized in 
Tabic |H 



First, we consider the general structure of the single- 
particle Hamiltonian in the two-band language diagonal 
in the electron operators, i.e. the general form of the ki- 
netic energy. This part of the Hamiltonian is understood 
in terms of hopping. 



1. Spin-independent terms 



The spin-independent part is given by 

^ = E E " kQQ ' C aks C Q'ks 



(2) 



a, a' k,s 



where for ordinary hopping the energy term Skaa' = 
Skaa^QQ' is diagonal in the band index. We may explic- 
itly write 



H = E[(^ tra -^ + 4 ntcr )4 s c lks 

k,s 



„intra 



H £k ter )4ks C 2ks 



(3) 



Here, e k 3tra = ej^+Q is an intra-sublattice term, i.e. rep- 
resents hopping between sites of the same sublattice 
type. Then, ejf tcr = — ej^q is correspondingly an inter- 
sublattice term due to hopping between sites of different 
sublattices. It is now useful to introduce Pauli matri- 
ces (t , f) for the orbital space with which the matrix 
element in Eq. ([3]) simplifies to 



n / intra ,,\_u _i_ _ 

"kaa' ^ \ e k — h L ) T aa' "+" e 



inter 3 



(4) 



A sublattice dependent chemical potential (J-a(b) = A 4 ± 
A/i would lead to an additional spin-independent term 



2. Spin- dependent terms 

We now turn to the spin-dependent part of the Hamil- 
tonian, 



(8) 



a. a' k s,s' 



This single-particle Hamiltonian will be written as a ten- 
sor product of a spin and a band part. In the following, 
summation over repeated indices is implicit. 

First, we consider terms which are based on intra- 
sublattice contributions, connecting only sites of the 
same sublattice. These can be written as 

Iw = 4° • ® T° aa , + .4 1 • B ss > ® ri a , , (9) 

involving intra- and inter-band terms, f£ and f^. Anal- 
ogously the inter-sublatticc part is given by 



kao 



9k ■ °ss 



9k ■ °ss 



(10) 



These terms are important, if time-reversal and inversion 
symmetry are violated. Time reversal T and inversion / 
operate as 



and 



77 k a = -f-k 



Ifk = f-k 



(11) 



(12) 



and analogously for g£ . 

For illustration, let us look at a few generic examples. 
Zeeman coupling of all spins to a uniform magnetic field 
Hq is implemented by f£ = <?/iB-ffo and, correspondingly, 
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a staggered field Hq (opposite for the electron spins on 

the two sublattices) is represented as = g^HQ (anal- 
ogous to the sublatticc dependent chemical potential), 
both being on-site-coupling (intra-sublattice) terms. Ac- 
cording to Eq. ((TTj) , they introduce a violation of time- 
reversal symmetry. Spin-dependent hopping terms con- 
necting the same or different sublattices can be written 
as 

4° = AL ntra and & = ^» t , (13) 

respectively. More important for our subsequent discus- 
sion are "staggered" spin-orbit coupling terms which cor- 
respond to 

/k=C k ntra and .g k 2 = Ck tGr , (14) 
for intra- and intcr-sublatticc hopping, respectively 

3. Symmetry considerations 

We consider now some symmetry aspects, whereby the 
sublattice structure again plays an important role. We 
introduce Q as the generating point group and denote by 
Q' the subgroup of operations respecting the sublatticc 
structure. All other operations in Q \ Q 1 interchange the 
two sublattices. As we consider ccntrosymmctric crys- 
tals, we request that the inversion is contained in Q. 
However, inversion may or may not be contained in Q' . 
In the former case, the symmetry center lies on one of the 
two sublattices and the operations in Q \ Q' have to be 
accompanied with a translation undoing the interchange 
of the sublattices. This means that the space group of 
these crystals does not contain Q as a subgroup and is 
therefore non-symmorphic. In the latter case, the cen- 
ter of inversion lies between the sublattices. This can 
lead to both, symmorphic and non-symmorphic crystal 
structures as we will see in the following. 

As noted above, the diagonal terms of the single- 
particle Hamiltonian have a tensor product structure, 
consisting of the momentum-dependent spin part (ekC°, 
/k • a and <7k ' <?) and the band part expressed by the 
T-matrices. Therefore, we may classify these terms by 
means of irreducible representations of G, as R(E)R' . The 
symmetry operations g S Q on the momentum-dependent 
spin part act as 

f£k = e - k and 5 / k = D+f D - k , (15) 

where D~ is the corresponding operation of element g on 
a vector and on a pscudo- vector. For the band part, 
it is easy to see that r° and r 3 do not change under such 
an interchange of the sublattices, such that they belong 
to the trivial irreducible representation Ai g of Q . On the 
other hand, terms with r 1 and r 2 change sign under the 
interchange of sublattices and belong to an irreducible 
representation V specific to Q and the sublattice struc- 
ture. Depending on whether inversion is an element of 
Q' or not, V will be an even or odd representation. 





intra-sublattice inter-sublattice 


IR 


intra-band 


r° r 3 


Aig 


inter-band 


r 2 


r' 



TABLE I: The different band dependencies possible for terms 
in the Hamiltonian of the systems under investigation here. 
While t° and r 3 always belong to the irreducible representa- 
tion Ai g , the irreducible representation V' of the other two 
Pauli matrices depends on the symmetry operations, that 
have to be combined with a sublattice interchange to map 
the crystal onto itself. 

For illustration, we consider two specific examples for 
a lattice with tetragonal symmetry with Q = D^h, which 
will be discussed in more detail below. The first ex- 
ample has a sublattice structure such that the A- and 
-B-sublattices form alternating layers along the z-axis, 
which yields Q = (0,0,7r/c) and the primitive lattice 
vector interconnecting two sublattice points is (0,0, c), 
see Fig. HJa). The center of inversion lies in the mid- 
dle between the two layers, e.g. at (0,0, c)/2 and in- 
terchanges the two sublattices (a and c being the lat- 
tice constants in-plane and out-of-plane, respectively). 
In this case, the subgroup leaving the sublattices invari- 
ant is Q' = C4 V and T' = A2 U . The second example is 
a sublattice structure within each layer with the primi- 
tive lattice vector (a, a, 0) connecting the two sublattices, 
leading to Q = (n/a, n/a,0), see Fig.QJb). The inversion 
center lies within the layer on a lattice point belonging to 
one of the two sublattices. The subgroup retaining the 
crystal structure is C^h and V = A 2g of D^. 

B. Off-diagonal single-particle terms 

We now introduce the superconducting order param- 
eter which on the mean-field level leads to off-diagonal 
terms to the single-particle Hamiltonian. These terms 
can be classified in a very analogous way as the diagonal 
terms. It is illustrative to discuss first the pair wavefunc- 
tion 

*W = (CaWsC a >-Ws'), (16) 

which combines two electrons characterized by spin and 
band configuration. Note that the pair wavefunction de- 
scribes zero-momentum pairs for a = a' while for a =/= a' 
the pairs possess momentum Q as can be seen from the 
definition of the single-particle operators in Eq. ([1} . 

In order to formulate the off-diagonal terms in the 
Hamiltonian, we introduce now the (mean-field) gap 
function A* s a , (k) and write 

n' MF = A ^'( k ) c ik S c l<- ks < + h - c - (i7) 

k 

We use the standard notation of a scalar gap function 
ip(k) for spin-singlet and the vector gap function rf(k) 
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r+ 


^0,1,3 (k) 


<fe(k) 




1 


xkykz " ykzkx 


Ai a 




xk x k z ~\~ ykyk z 


Big 


kl-kl 


xkykz ~\~ ykxkz 


B 2 g 


k>x ky 


xkxkz ykykz 


Eg 


\kxk z 5 ky k z } 


\ikxk z j zkykz\ 


r~ 


V> 2 (k) 


^0,1,3 (k) 


Ai u 


- 




A 2u 






Biu 


kykz 


%>k x yky 


B 2u 


kz{k x ky) 


xky -J- y/^a; 


E u 


\k x 5 fcy } 


\ik x , zky\ 



TABLE II: Basis functions belonging to the different irre- 
ducible representations of D4h with SOC for the different 
gaps. 

for spin-triplet pairing. The gap function has to satisfy 
the Pauli principle to change sign under exchange of the 
two electrons: 

A»£(k) = -A&(-k). (18) 

For a single-band superconductor this requires that 
V>(-k) = tp(k) and d(-k) = -d(k). 
We express the gap function as 

A£,(k,a) - hU%° + d a (k) • c] ss , ® r° a „ (19) 

where we define <r° = icr 37 and <f = ioo v . For intra- 
sublattice Cooper pairing originating from interactions 
between electrons on the same sublattice, intra-band cor- 
responds to a — and inter-band pairing to a — 1. Anal- 
ogously, inter-sublattice pairing for intra-band pairs takes 
the index a = 3 and for inter-band pairs a = 2. Note that 
for a = 0, 1 and 3 the scalar (vector) gap function is an 
even (odd) function of k, while it is opposite for a = 2, 
as required by Eq. (fl8|). The case of a = 2 is special 
in the sense that the band part of the pairing state is 
antisymmetric under exchange allowing both momentum 
and spin part to be simultaneously symmetric or anti- 
symmetric. 

As in the case of the diagonal part we can classify the 
symmetry for the tensor product characterizing the pair- 
ing state (gap function). Thus, we consider again the ir- 
reducible representations R s <g> R' s of the generating point 
group Q. The representations R' s correspond again to the 
ones of the r matrices as given in Tab|TJ The represen- 
tation R s is based on the internal (spin and momentum) 
structure of the Cooper pair, given in TabHTlfor the case 
Q = Dih j which we will use in the following. 

Together with the representations of the diagonal part 
it is possible to see which types of pairing states would 
couple together due to the staggered form of spin-orbit 
coupling. Such couplings occur due to combination of 
the diagonal and off-diagonal part in the linearized gap 



equations as a form of selection rules^. Thus, having the 
pairing state R s ® R' s also a state of symmetry R s ® R' s 
would be coupled through a diagonal term of symmetry 
i? £§> R' , if it appears in the decomposition of 

(R <g> R') x (R s ® R' s ) = (Rx R s ) <g> {R' x R' s ). (20) 

This allows to classify all possible interdependent pairing 
states within a given crystal lattice symmetry. 

In the following, we will analyze three different tetrag- 
onal crystal lattices with generating point group D^h and 
elaborate on the way of analyzing the influence of stag- 
gered types of spin-orbit coupling due to local inversion- 
symmetry breaking on superconductivity from a symme- 
try point of view. 



III. STACK OF INVERSION SYMMETRY 
LACKING LAYERS 

Our first example is a tetragonal crystal lattice, whose 
staggered form originates from a sublattice structure of 
alternating layers. The basic unit, the layer, violates in- 
version symmetry by the absence of reflection symmetry 
z —z (z: the four-fold rotation axis of the tetragonal 
crystal). This type of non-centrosymmctricity yields a 
Rashba-type spin-orbit coupling in each layer, 

H soc = J2(^-^M s c ks , (21) 

k 

with Ak = a(x sin k y — ysmk x ). The sign of the Rashba 
coupling a is opposite for the two sublattices, i.e. alter- 
nates from layer to layer (see Fig. HJa)). The two bands 
resulting from this feature are related by Q = (0, 0, n) 
taking from now on all lattice constants to unity. 



A. Symmetry considerations 

The crystal lattice has the tetragonal D^h point group 
with full inversion symmetry, taking the center at a sym- 
metry point between the layers. The elements of D^h 
are divided into those transforming within the layers and 
those interchanging the sublattice: 

G intra = { E ,2C 4 , C 2 ,2a v ,2a d } = C Av , (22) 
G intcr = {2C 2 ,2C^,I,a h ,2S i }, (23) 

using the standard notation of Ref. 0- From this we con- 
clude that R' — Aiu, which is the one-dimensional irre- 
ducible representation with +1 for all elements of G mtra 
and —1 for all elements of G mtcr . Considering now the di- 
agonal single-particle part of the Hamiltonian, we find for 
the spin-independent hopping terms the standard repre- 
sentations A\ g (£> A\ g . On the other hand, the staggered 



5 



spin-orbit part consists only of the intra-sublattice (in- 
plane) Rashba-like coupling for which transforms ac- 
cording to A 2u and corresponds to = Ck ntra of Eq. (fT4|) . 
This leads to the representation A 2u <8> A 2u . 

For this system, spin-orbit coupling mixes pairing 
states in the following way, 

{A 2u ®A 2u )x(R s ®R a <) = (A 2u xR s )®(A 2u xR' s ), (24) 

i.e. states of opposite parity can be mixed, as is generally 
the case in non-centrosymmetric systems. Note that this 
also implies that intra-band pairs mix with inter-band 
pairs. 

Looking first at intra-sublattice (intra-layer) pairing 
states, we consider the example of the (even-parity) s- 
wave spin-singlet state, which has for intra-band pairing 
the representation Ai g ®Ai g while it belongs to Ai g ®A 2u 
for inter-band pairing. The mixing occurs as follows, 



A 



is 



Alg O A 2 



A, 



A lg <g>A 2u «-» A 2u ® A lg , 



(25) 



whereby the admixed states have always opposite parity 
(A 2u ). Using TablTIlwc write the two types of states with 
intra-sublattice pairing as 

A(k) =i) <; ®T + d 1 (k y <; x -k x <;y)®T 1 , . . 
A(k) =Vi?° &T 1 + do{k y <; x -k x $y)®T°, 

which mix the spin-singlet and triplet configurations. 
Note that the same scheme also applies for other pairing 
states, e.g., a d-wave state beloning to Bi g (£>Ai g couples 
to a spin-triplet pairing state belonging to B 2u (g> A 2u . 

Next we consider inter-sublattice (inter-layer) pairs, 
starting with s-wave intra-band states, corresponding 
again to A\ g <E> Ai g with the admixed A 2u <£> A 2u . On the 
other hand, the inter-band (even-parity) "s-wave" state 
(Ai g <g> A 2u ) has a spin-triplet configuration and couples 



is 



to the intra-band odd-parity spin-triplet state A 2u 
as the t 2 matrix is involved (Tab]]}. 

For the two possible inter-sublattice pairing states we 
find the gap functions 



A(k) = 



(27) 



on a tight-binding band structure, whereby each layer is 
considered as a simple square lattice. We use the two- 
band formulation and write the single-particle part of the 
Hamiltonian as 



k 



ks' 



(29) 



with electron operators as defined in Eq. (jTJ) with Q = 
(0,0,7r). Intra-layer hopping is taken into account be- 
tween nearest and next-nearest neighbors and inter-layer 
only between nearest neighbors, which leads to 



-intra 



-inter 



-2t xy (cos k x + cos k y ) — it' xy cos k x cos k y 



-2t? cos ky 



(30) 

These contribute to the spin-independent part. The spin- 
dependent part originates from the staggered Rashba- 
typc spin-orbit coupling, which we take only in the 
nearest-neighbor form as in Eq. (|2ip yielding 



/k = a (% sm ky ~ V sm kx) 



(31) 



following Eqs. © and ([TP]). 

It is convenient for the following to use the formula- 
tion by means of Green's functions, which for the non- 
interacting case can straightforwardly be calculated by 
inverting the (4 x 4) matrix (iuj n a° s , ®r° Q , - Hk aa '0-° s , - 

yss' \ 
kact'/' 

G (k,w„) = G Q+ (k,oj n )<7° <S)t° 

+ G -(k, u> n )(f k • <g> r 1 + i w <j° g> r 3 ), (32) 



where 



G ±(k,uj n ) = \(- 1 —± 



2Viw n -^ +: k iu n - £_ 



/k=/;vv / i4 l i 2 +(4 nter ) 2 



(33) 



(34) 



and 



and 



A(k) = d 2 (kyk Z ^ ~ k X k Z ^) ® T 2 



(28) 



which remain in either the spin-singlet or spin-triplet 
channel. Due to the sublattice structure, however, al- 
ways inter- and intra-band states are mixed. 



B. Microscopic consideration 

To illustrate the symmetry-based aspects from a mi- 
croscopic point of view we introduce here a model based 



ek=e k ltc VVI/ k 1 l 2 + ( £ k tcr ) 2 - 



(35) 



In Eq. (|33[) . the two (spin- independent) band energies are 
given by 



£±,k S = £±, k = 4 ntra - m ± V 1/kT + (4 nter ) 2 - (36) 

We now turn to the problem of superconductivity by 
introducing a pairing interaction of the general form, 

^ = T7 VaP^u 4 (^' ^ ) C t*ks C l 9-ks' C ;i-k's3 C Vk's4 - 



k,k' 



(37) 
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We parametrize the matrix element in the notation used 
for the single-particle terms, 



c s r 4 (k,k') 



.,(<*) r,/,( a ) 



x [itilWC^W (38) 



where V'mn(k) have the symmetry of the gap functions 
tabulated in Tab|n] For a more detailed analysis of the 
structure of such an interaction see appendix |XJ This 
pairing interaction incorporates both, coupling of the 
intra- and inter-sublattice type. For simplicity, we re- 
strict ourselves to interactions including only nearest- 
neighbor coupling in the real lattice. This limits the 
classification of pairing states as can be seen in Tab lllll 
compared to the more general Tabim 

With the Hamiltonian and the non-interacting Green's 
function introduced above it is possible to analyze the 
superconducting instabilities in detail by resorting to the 
standard framework of the Gor'kov equations^. The lin- 
earized gap equation reads 

a$(io = -t^T^T, E ^!i S4 ( k > k ') 

f_L , V LO n k' S3 } S4 

x [G (k', w „)A(k')G^(-k',- w „)]^, (39) 

where all the Green's functions as well as the order pa- 
rameter are 4x4 matrices. This gap equation is analyzed 
in the following for the two cases of a leading instability in 
the intra-sublattice and the inter-sublattice pairing chan- 
nel, respectively. 



1. Intra-layer interaction 

We use the nearest-neighbor interactions derived in ap- 
pendix for the intra-sublattice case, which, following 
Eq. (|26|) . lead to gap functions of the form 



A(k) 



^o(k)<? ®T° + d 1 (k)-<r®T 1 

V-iik^^+doikj-^T 



for which we insert from Tab lllll 



V>n(k) = 1p n (c<38k x +COSfcj,), 

d„ (k) = d„ (x sin k y — y sin k x ) 



(40) 



(41) 
(42) 



with n = 0, 1. It is easy to see that the gap functions (|4"0)) 
couple within the linearized gap equation (|39[) indeed in 
the way anticipated above, using the intra-layer interac- 
tion given in the appendix in Eqs. (|A19[) and (I A20I) : 

ifo(k) = -T^«+ k ,{[Go + G 0+ + Go-Go_]Vo(k / ) 

n,k' 

+ [Go+Go- + G -G +]/ k < • di(k0}, (43) 





intra-sublattice 


inter-sublattice 




^0,1 (k) 


ib 3 (k) 




1 pnQ tc -4— rr»Q tc 
± , 1 ' ji 1 nix Lua ^y 


cos k z 


Ri 


pn<* I? — pns 1? 

v ' ro rn x ^ \ > . 1 






do,i(k) 


^ 2 (k),d 3 (k) 


Aiu 


x sin k x + y sin k y 


1 sin fc z 


A 2u 


x sin k y — y sin k x 


sin k z 


Bin 


x sin k x — y sin k y 




B2u 


x sin k y + y sin k x 




E u 


{z sin k x , z sin k y } 





TABLE III: List of different basis functions for a crystal struc- 
ture with an alternating stack of mirror-symmetry lacking 
layers that are supported by nearest-neighbor intra-sublattice 
and inter-sublattice interaction, respectively. Note that with 
the restriction that only nearest-neighbor pairing is consid- 
ered the "inter-sublattice" case does not include any even- 
parity spin-triplet d2(k)-states. 

di(k) = -T^ l)kk ,{[G 0+ G 0+ + G -G -]di(k') 

n,k' 

+ 2Go-G -{.M/k< • di(k')] - di(k')} 

+ [Go+Go- + Go_G 0+ ] jWo(k')} (44) 
and, analogously, 

d (k) = -T^« kk ,{[G 0+ G 0+ + Go- Go-] d (k') 

n.k' 

+ 2G -Go-{.M/k< • do(k')] - do(k')} 

+ [Go+Go- + G -Go+]/ k 'Vi(k')}, (45) 

lh(k) = -T^« kk ,{[G +G + +G -G -]Vi(k') 

-2(e k 2 Go-G -Vi(k') 
+ [Go+Go- + Go-Go+]/ k ' ■ do(k')}. (46) 

Here, we have introduced the short notation Go± = 
Go±(k, u n ) and Go± = Go±(— k, — w n ). It is also obvious 
now that it is the spin-orbit coupling term, represented 
here through / k , which yields the coupling between even- 
and odd-parity pairing states. 

2. Inter-layer interaction 

Turning to the inter-sublattice (inter-layer) pairing the 
situation becomes more intricate due to pairing in the 
anti-symmetric band channel (r 2 ). We write the gap 
function as 



1 ' I V 3 (kK 3 ® r 3 + d 2 (k) ■ f® r 2 



(47) 
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III 


dk = yk x - xky — 
_ dk = §k x + xky - - - 
dk = xk x + yky 
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FIG. 2: Suppression of the transition temperature of the 
intra-band gaps due to the antisymmetric SOC. With in- 
creasing inter-layer coupling, the suppression is weakened, 
t z — 0,0.1t,0.2t from bottom to top. 

where for nearest-neighbor pairing only we may use 

$j(k) = ^ 2 sinfc z , (48) 
V> 3 (k) = ipscoskz, (49) 
d3(k) = d-ii sin k z , (50) 

while d 2 (k) = 0, following TabHID To write the lin- 
earized gap equation we use the inter-layer pairing inter- 
actions (|A21|) and (|A22j) to find 

lfe(k) = -T^+ k , {[G 0+ G 0+ + G _G _]^ 3 (k') 

n.k' 

-2G -Go-/^V3(k')} (51) 

and, in the same way, 

d 3 (k) = -T^ ? ; kk ,{[G 0+ G 0+ +G _G _]cr 3 (k') 

n,k ; 

-2G -Go_[/ k -d3(k')]/k'}. (52) 

Note that for the gap equation there is no mixing within 
this approximation. 

C. Discussion 

We first consider the intra-layer pairing channels look- 
ing at Eqs. (|43|) - (|4T))) . Restricting our analysis to the 
terms diagonal in each gap function only, we find that 
even-parity intra-band pairing is essentially unaffected 
by spin-orbit coupling. More interesting are the intra- 
band odd-parity gaps which suffer suppression unless the 
d vector is parallel to /k. This property is known from 
non-centrosymmetric superconductors^ and should be 
fully transferable to the case of completely decoupled lay- 
ers. Our analysis shows that including inter-band pair- 
ing combined with inter-layer hopping [t z ^ 0) reduces 



the pair breaking effect of spin-orbit coupling, as can be 
seen in Fig. [2l There, we plot T c versus the spin-orbit 
coupling strength a for three different values of the inter- 
layer hopping values. Naturally, the case of do(k) || /k is 
unchanged (solid line in Fig. [5]). 

Considering the inter-band gaps [^i(k) and di(k)], 
which correspond to finite-momentum pairing, it is easy 
to see that they have the same gap equations as the 
intra-band gaps, if we decouple the layers (t z = with 
£Tk = 0) and set the spin-orbit coupling to zero (a = 
with /k = 0). The reason is that the finite momen- 
tum Q = (0, 0, 7r) yields an alternating phase of and n 
from layer to layer, which is irrelevant for decoupled lay- 
ers. Therefore finite-momentum pairing is suppressed by 
inter-layer hopping, which introduces the disadvantage 
of an inter-layer phase shift to the energy balance, as we 
can see in the plot of T c in Fig. [3l On the other hand, 
adding spin-orbit coupling helps the intra-band pairing 
to slightly recover T c . 

The trends discussed so far show that a strong inter- 
layer coupling moves the system further away from the 
parity-mixing (spin singlet-triplet mixing) as compared 
to a real non-centrosymmetric superconductor. Inter- 
layer coupling "recovers the inversion symmetry" gradu- 
ally. 

A further interesting aspect occurs, if the inter-layer 
pairing is the leading instability. While the system has 
the full inversion symmetry in this case, the spin-orbit 
coupling acts as pair-breaking for the spin-singlet channel 
(see Eq. (fSTj) '). On the other hand, spin-orbit coupling 
does not affect the spin-triplet gap (^(k), since ^(k) _L 
/k for all k. Thus, in this case the spin-orbit coupling 
can be important in influencing the pairing symmetry 
in favor of a spin-triplet state. Note, however, that the 
structure of the pairing interaction remains the major 
deciding element for the pairing symmetry. 



IV. INVERSION SYMMETRY LACKING 
BONDS 

The second example studied here is motivated by the 
layered pcrovskite crystal structure with tetragonal sym- 
metry, known for some transition metal oxides. The sub- 
units are oxygen octahedra, where six oxygen ions enclose 
a transition metal ion. In some systems, these octahedra 
rotate around the crystalline z axis leading to a stag- 
gered pattern of rotation, i.e. neighboring octahedra in 
the xy plane rotate in opposite direction [see Fig. QJb)]. 
Such features are known in the bilayer Sr3Ru2 07 and 
in Sr2Rh04 or Sr2lr04, to mention a few examples.—^ 
This lattice distortion shifts the in-plane bond oxygens 
to off-center positions, and thus, leads to a breaking of 
inversion symmetry on each of these bonds. Again, we 
arrive at a form of staggered spin-orbit coupling fitting 
into the scheme developed above. Since inversion sym- 
metry within the sublattice is retained in this structure, 
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FIG. 3: Change in the transition temperature of the inter- 
band gap with di (k) = x sin k y — y sin k x as a function of 
the inter-layer hopping t z for different spin-orbit coupling 
strengths. 



there is no even-odd mixing at all for this type of crystal 
structure. However, the spin-orbit coupling has an in- 
fluence on the direction of the d-vector for a spin-triplet 
pairing state. Before this is studied in detail with the 
help of the linearized gap equation, we again start with 
a symmetry analysis. 



A. Analysis of symmetry 

Without the rotation of the 06 octahedra, the point 
group of a single layer is D 4 h as in the above example. 
With the rotation, the unit cell is doubled in plane as 
depicted by the light and dark lattice sites in Fig. QJb) . 
The corresponding Q-vcctor is Q = [tv,tx, 0). 

We separate again the symmetry operations within 
D^h which turn each sublattice into itself (G mtra ) while 
the remaining operations (G lnter ) exchange the sublat- 
tices, 



C?intra — {E, 2C±, C%,I, 2Si, (J h} — C±h, 

G intor = {2C 2 ,2C 2 ',2a v ,2a d }. 



(53) 



The representation of D^h changing sign for all elements 
of G intor is T' = A 2g . 

Analogous to the previous example, the terms in the 
Hamiltonian can be characterized with respect to their 
behavior under sublattice interchange. According to the 
above symmetry analysis, the terms that change sign be- 
long to the irreducible representation V = A 2g for this 
structure and a symmetry-reducing term in the Hamil- 
tonian has to be of A 2g ® A 2g symmetry. The staggered 
spin-orbit coupling derived in Ref. [ll] is of this symmetry 
with 



<7 k = 2a (cos k x + cos k y )z = C k n1 



(54) 



It is again possible to find the symmetry allowed cou- 
plings between different gap functions by factorizing the 
gap function in terms of spin and orbital degrees of free- 
dom to find the general relation for mixing pairing chan- 
nels, 



V® A Xg <-> (r x A 2g )®A 2gi 
which, for instance, leads to 



.4 



19 



B 



lg,u ' 



E, 



9," 



A Xg o A 2g ®A 2g , 
Ai g <-> B 2g ^ u ®A 2g , 
Ai g <-» E gtU <g> A 2g . 



(55) 



(56) 
(57) 
(58) 



As mentioned above, spin-orbit coupling here does not 
mix states of different parity. 



B. Analysis of instability 

A better understanding of the consequence of symme- 
try properties and of the influence of the spin-orbit cou- 
pling on the different superconducting states can be ob- 
tained by analyzing the linearized self-consistency equa- 
tion for the gap (|39[) . The non-interacting Hamiltonian 
is the same as in Section Mil with the only difference 
that the spin-dependent term here uses g£ as given in 
Eq. (|54|) . Note that we restrict ourselves to the single- 
band case and ignore the aspect of degenerate d orbitals 
of transition metal ions in the examples mentioned above. 

The non-interacting Green's function is given by 

G (k,w n ) = G Q+ (k,u n )<j° ®t° 

-Go-(k,w n )(g k CT z (g>T 2 -e k CT°®T 3 ), (59) 



with 



G ±(k,u) n ) = -( 



± 



1 



and 



k.s 



.9k 



2 iuj n - £ +:k iu n - £_ 

(5kWv / lflk| 2 + (4 nter ) 2 : 

^V^T + (4 nter ) 2 , 



- M±v /|ff k 2 l 2 + (£k 



~inter\2 



(61) 



(62) 



(63) 



and has even parity. 



Note that we again distinguish hoppings connecting dif- 
ferent sublattices (inter-sublattice, e k ator ) and the same 
sublattice (intra-sublattice ejf tra ) with the former in- 
cluding nearest-neighbor and the latter including next- 
nearest-neighbor hopping. In the following, the two cases 
of an intra-sublattice and inter-sublattice pairing inter- 
action arc again discussed separately. 



1. Intra- sublattice pairing 

For a leading interaction of intra-sublattice type, the 
gap is analogous to the form given in Eq. (|40[) . As men- 
tioned, the intra- and inter-band gap functions have the 
same parity and spin configuration. The case of even- 
parity pairing can be illustrated with the example of the 
intra-band pairing state 



Vo(k) = 
of D 4h . 



sin fc T sin fc„ 



(64) 



belonging to Bi g of D^h- The corresponding admixed 
state according to Eq. ([57)1 is the i?i g -state 



ipi (k) = cos 2k x — cos 2k y 



(65) 



Note that i/jq (k) here is based on next- nearest- neighbor 
pairing, while ipi (k) originates from an interaction on 
sites separated by lattice vectors (2a, 0). 

For the odd-parity channel it turns out that the x and 
y components of the d vector for intra- and inter-band 
pairing states mix in the form 



dg(k) 

dm 



-d?(k) \ 
dm ) 



which is a result of the decomposition of T < 
leads to combinations of gap functions like 



(66) 



This 



A\ u : do(k) = xsm(k x + k y ) + ysin(k x - k y ) «-» 

A 2u : rfi(k) = xsin(k x - k y ) - ysm(k x + k y ), (67) 

representing one example of a pairing state classified 
within the representation of T)^h and arising from next- 
nearest-neighbor interaction. The z component of the d 
vector is conserved in the mixing of inter- and intra-band 
pairing. This is fully compatible with the classification of 
spin-triplet pairing states in a tetragonal crystal lattice. 
The state belonging to E u then yields for next-nearest- 
neighbor pairing, 

d (k) = {i sm(k x + k y ),z sm{k x - k y )} (68) 

mixing with 

d\ (k) = {— z sm(k x — k y ), z sm(k x + k y )}, (69) 
which also lies in the representation E u . 

2. Inter- sublattice pairing 

For inter-sublatticc pairing interactions, we again con- 
sider a gap of the form given in Eq. (|47j) . If we only con- 
sider pairing within the xy plane, the spin-singlet pairing 
states "03 (k) only appear in the one-dimensional repre- 
sentations of D^h, while the states ^(k) are in the two- 
dimensional representation E u (the respective others re- 
quire that the gap function changes sign under the opera- 
tion z —> —z). Therefore, following Table ITVl we find that 



r+ 


Vs(k) 


<fe(k) 


Alg 


cos k x + cos k y 


- 


Mg 




z (cos k x + cos k y ) 


Big 


cos k x — cos k y 




B 2 g 




z (cos k x — cos k y ) 


Eg 




{£(cos k x + cos ky), y(cos k x + cos k y )} 


T~ 


V> 2 (k) 


d 3 (k) 


Ai u 


- 


x sin k x + y sin k y 


A 2u 




x sin k y — y sin k x 


Bi u 




x sin k x — y sin k y 


B 2n 




x sin k y + y sin k x 


E u 


{sin k x , sin k y } 


{z sin k x , z sin k y } 



TABLE IV: Basis functions belonging to the different irre- 
ducible representations of D4h supported by in-plane nearest- 
neighbor interactions. 



corresponding spin-triplet components d(k) _L z remain 
independent. Only the E g ^ u spin-triplet state (d(k) || z) 
mixes with the spin-singlet states. 

We consider first the case d(k) _L z yielding the follow- 
ing linearized gap equation, 

#»(k) = -T^4« kk ,(G 0+ G 0+ + G -G _)d^(k') 

n.k' 



and 

d x ' v 



(70) 



2 (k) 



-T^4 v + k ,(Go+G + - Go-Go-^fk'), 



(71) 

where we used again the short-hand notation Go± = 
Go±(k',w„) and Go± = Go±(— k', — uj n ) and u kk , are de- 
fined in App. [S] 

On the other hand, the z component mixes with a 
scalar gap function, 




= -r£4« kk , [M(k')]( 

n,k' V 



df(k') 

MX) 



(72) 



for the odd-parity and, similarly, for even-parity gap 
functions, 




n,k' V 



^ 3 (k') 
df(k') 



(73) 



The matrix in Eqs. (|T2"|) and ([73")) is given by 

Mn(k) = Go+G + +G0-G0- -2^G -G -, (74) 
Af22(k) = G0+G0+ + Go-Go- — 2e k Go_Go-, (75) 
M X2 (k) = 2i5 k £ k G -Go_ = M| x (k). (76) 

Performing the sums over the Matsubara frequencies, 
we first discuss the uncoupled x and y component of 
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the d vector and choose the odd-parity gap functions, 
c?3(k) = (A^.t + A y _y)smk x obtained for a nearest- 
neighbor interaction [a degenerate solution of the lin- 
earized gap equation is d 3 (k) = (A^x + A v __y) sin k y }. 
Eq. ([70)) yields then the standard BCS equation deter- 
mining T c which is degenerate for both x and y compo- 
nent, 



with 



l = -y££^tanl/^ k 



k' a=± 



2T 



(77) 



We turn now to the even-parity gap function ^(k) 
for which we assume an extended- s- wave form, ^(k) = 
(A x + x + A y + y)(cosk x -\- cos fcy), a result of the nearest- 
neighbor interaction. This leads to the equation for T c , 



^ ^ (cos k' + cos k' ) / ' E a v> ! 



2T 



originating from Eq. (|71[) . 

These equations should be compared with the corre- 
sponding equations ([72)1 and (|75|) . For the odd-parity 
case with the nearest-neighbor coupling approach and the 
gap functions d\ (k) = Af. sin k x and ip2 (k) = Af_ sin k x , 
we obtain from Eq. (1721 



A? 
A: 



L +L\ 
—iL 3 L Q 



Li 




(79) 



Summing again over the Matsubara frequencies, we can 
express these matrix elements as 



Lo = 


-yj^sin 2 

k 


k x Si{\s), 




Li = 


-y^sin 2 

k 


k x gl[S 2 (k) - 


Si(k)], 


L2 = 


-y^sin 2 

k 


k x el[S 2 (k) - 


Si(k)}, 


Ls = 


-yj^sin 2 


^ffkEk^Ck) 


-5 x (k)] (80) 



with 



Sl (k) = y _L t anh(^) 
t± 2C a , k \2T ) 

Sa(k) = J2 



(81) 



^ 2(£ mtra _ ^ 



The analogous result can be obtained for the even-parity 
case with df(k) = A^(cosfc :E + cosfc y ) and "03(k) = 
A^cosfcz + cosfc y ), leading to 



A3 




(83) 



L o = 


— V >(COS K x - 
k 


hCOS fcy^ 




Lt = 


— y^^(cos fc x - 

k 


hcos k y 


) 2 .g 2 [5 2 (k)-5!(k)], 


Lt = 


~y^(cOS 
k 


hCOS fcj, 


| a ^[5 a (k)-5i(k)], 


Lt = 


— yy^(cos k x - 


hcos k y 


) 2 .9k£k[ft(k)-5i(k)].(84) 



The instability condition for Eq. ([79"f and (1531 are given 
by the eigenvalues 



A s ± = LI + 1{L\ + L s 2 ) ± i v /(Lf-L5)2 + (2L|)2 (85) 

reaching AJ. = 1 for both even and odd parity with s = + 
[Eq. ([75)1 ] and s = - [Eq. ([83)l]. Note that the instability 
condition for the a; and ?/ components in Eq. (|77|) and 
([78)1 correspond to 



(86) 
(87) 



Ai = Lo = 1 
for odd-parity pairing [Eq. (|77[) ] and 

A 2 = L+ + L+ + L+ = 1 



for even-parity pairing [Eq. (|78[) ]. It can be demonstrated 
easily that Lf,L| < and L\L 2 > (L|) 2 . We now use 
the resulting inequality, 

(^-L^) 2 + (2^) 2 <(^+LJ) 2 (88) 

and A:j_ > Ai to obtain the relation: 

L s > \ & + > X S L > L & + L{ +L S 2 . (89) 

From these relations we are able to show for the odd- 
parity states 



A x — L > A ± 



(90) 



such that the instability leads to a state with the d vector 
perpendicular to the z axis as described by Eq. |70)) . On 
the other hand, for even-parity pairing the inequality, 



A± > A 2 — Lq + L\ + L'2 



(91) 



favors the state with spin-singlet and -triplet mixing 
where the d vector points in z direction as described by 
Eq. J83)). 



V. STAGGERED NON-CENTROSYMMETRIC 
PLAQUETTE STRUCTURES 

As a further application we turn to a system with two 
sublatticcs each lacking inversion symmetry, motivated 
by the crystal structure of some of the iron-pnictide su- 
perconductors. There, a single FeAs layer consists of Fe 
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(b) 



r. 



FIG. 4: (a) Top view of the basic FeAs crystal structure. The 
open circles denote As ions lying below the plane while the 
crosses denote ions above the plane, (b) Top view of one of 
the sublattices rotated by 45 degrees for an easier analysis of 
the hopping Hamiltonian. 



ions forming a square lattice with As ions sitting in ev- 
ery center of the squares. As is depicted in Fig.@Ja), the 
As are shifted out of the Fe plane in a way as to built 
distorted tetrahedral cages around the Fe sites. Due to 
the arrangement of the As sites, this structure can again 
be described with two sublattices of checker-board type. 



A. Analysis of symmetry 

As in the first example of inversion-symmetry-lacking 
layers, this crystal possesses a center of inversion that is 
located between the sublattices. Here, however, the crys- 
tal structure is non-symmorphic. Taking the symmetry 
center on one of the Fe sites, we can again separate the 
symmetry operations within D^h leaving the sublatticc 
structure invariant and the ones interchanging the sub- 
lattices, 



G intra = j C 2 ,2S i; 2<J d } = D 2d , 

G intor = {2C 4 ,2^',/, a h ,2a v ,}. 



(92) 



As G lnter contains inversion, the representation of D±h 
that changes sign for all elements of G mter is again odd, 
namely V = B\ u . 

Following the characterization of the terms in the 
Hamiltonian introduced above, a symmetry-reducing 
term has to be of B\ u £§) B\ u symmetry. Such a term 
is microscopically derived in appendix [B] for a simplified 
orbital structure on the Fe sites and has the same or- 
bital and spin structure as the one in the first example, 
however with 



/ k = a(x sin k x cos k y — y cos k x sin k y ). 



(93) 



The symmetry allowed couplings between pairing states 
of different symmetry can again be found factorizing the 
gap functions according to 



leading, for example, to 

(Aig <g> Ai g ) o {B lu ®B lu ) (95) 
{A lg <E>B lu ) o (B lu ®A lg ). (96) 

This yields again states of mixed parity. 



B. Microscopic considerations 

The Hamiltonian describing this system has the same 
structure as the one encountered in Sec. IIII[ but with 
different dispersions given by 



£ mtra = -fa' C OS k x COS ky , 

4 ntcr = -2t{cosk x + cos ky), 



(97) 
(98) 



r®ii 9 n(rx b 1u ) <g> B lqj 



(94) 



and /k as defined in Eq. (|93[) . The linearized gap equa- 
tions have therefore the form of Eqs. (|4"3"| - (|4"6"|) for the 
intra-sublattice pairing and Eqs. (fS"Tj) and J52]) for the 
inter-sublattice pairing. Note that the difference in the 
crystal structure has more drastic consequences, since 
the possible pairing terms in the interaction for intra- 
and inter-sublattice interactions allow now for different 
gap functions as summarized in Table fVl Obviously, the 
leading pairing channel has to be at least a next-nearest- 
neighbor interaction to allow for spin-singlet to spin- 
triplet mixing for the intra-sublattice pairing. Therefore, 
we find for example a coupling of the two gaps 

Ai g : ipo (k) = ipo cos k x cos k y -h- 

B± u : di(k) = di(j)svak x cosk y — xshik y cosk x ). (99) 

Considering inter-sublattice pairing, the dominant 
channel is modeled by a nearest-neighbor interaction. 
Regarding first spin-triplet pairing, spin-orbit coupling 
lifts the degeneracy of the spin configuration. Obviously, 
states belonging to the representation E u with d || z 
are according to Eq. (f52"| unaffected and would yield 
the highest transition temperature. Any state in the 
other representations (A\ u , A 2u , B\ u , B 2u ) would have a 
reduced T c . On the other hand, the spin singlet pairing 
channels based on inter-sublattice interactions are gener- 
ally suppressed irrespective of the representation. 



VI. CONCLUSION 

For crystal lattices, where inversion symmetry is bro- 
ken in a regular, but non-uniform (unit-cell multiplying) 
pattern, the multi-band structure of the reduced Bril- 
louin zone renders the classification of the superconduct- 
ing order parameter in terms of standard spin-singlet 
and spin-triplet insufficient. On the level of the normal 
state electronic properties this is imprinted by spin-orbit 
coupling whose structure is closely connected to lattice- 
symmetry details. For such systems, the usual connection 
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intra-sublattice 


inter-sublattice 




1, cos k x cos k y 




cos k x + cos fc H 


Big 


- 




cos k x — cos fc H 


B 2g 


sin k x sin k v 




- 


Aiu 


x sin k x cos k y + y sin k y 


cos k x 


x sin k x + y sin fc H 


A 2u 


y sin k x cos k y — x sin k y 


cos fca; 


£ sin fc y — y sin k x 


Bin 


x sin k x cos k y — y sin fc y 


cos k x 


x sin k x — y sin fcy 


B 2 u 


y sin fca; cos k y + x sin fc y 


COS fca: 


i sin fc y + y sin fc^ 


E u 


{2 sin k x cos fc y , 2 sin k y cos fc^ } 


{2 sin k x , 2 sin 



TABLE V: Lowest order basis functions supported by intra- 
and inter-sublattice interactions on the lattice considered in 
Sec. [V] i.e., for on-site and nearest-neighbor interactions, re- 
spectively. In order to allow for a spin-singlet to spin-triplet 
coupling, an interaction between next-to-nearest neighbors 
has to be considered. 



between even (odd) parity in momentum space and spin- 
singlet (spin-triplet) configuration is lost in many cases, 
although the overall system is centrosymmetric. For a 
full classification in terms of the crystal symmetry, also 
the structure of the gap in band space has to be taken 
into account. 

In this paper, we have studied two main classes of a 
local lack of inversion symmetry with a two-sublatticc 
structure, whereby the complete lattice possesses a center 
of inversion. This corresponds to a doubling of the unit 
cell leading to a two-band description. The first class of 
lattices is characterized by the property that each sub- 
lattice has broken inversion symmetry. This manifests 
itself in the symmetry group G" ntra which does not in- 
clude the element of inversion /. This is the case for our 
first (Sect lIIl[ ) and last (Sect|VJ example. The situation 
yields singlet-triplet mixing which is characterized by the 
lattice specific representation V having odd parity and 
occurs in the connection with intra-sublattice pairing. 

In the other case, the sublattices retain separately in- 
version symmetry, i.e. G lntra contains /, while the links 
connecting the sublattices lack inversion symmetry. The 
corresponding representation T' has even parity, which 
also determines the structure of the spin-orbit coupling. 
In this system, it is inter-sublattice pairing which mixes 
spin-singlet and spin-triplet pairing while the parity re- 
mains fixed. Here obviously spin-configuration and par- 
ity are not anymore tied together. 

This new classification scheme can be important to de- 
termine which pairing states can be stabilized. This can 
be particularly useful, if questions concerning the degen- 
eracy in spin space have to be answered. The new states 
and electronic structures may have an impact on the way 
superconductors couple (Josephson effect) and how the 
superconducting state reacts on external magnetic fields. 
These topics will be discussed elsewhere. 
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Appendix A: Structure of Interaction 

In this appendix, the structure of a general density- 
density interaction in a crystal with a two-site unit cell 
is analyzed. The generalization to other types of inter- 
actions, e.g., a spin-spin interaction, is straightforward. 
Our starting point is a real-space formulation of the in- 
teraction, 

n' = E E Vijn is n js , = E E ^'iSv V c i* ,:Ali 

i,j s,s' i.j s,s' 

with Vij the interaction strength between the lattice sites 
i and j. Note that for the special case of i = j (on-site 
interaction) the spin sum only runs over s =/= s' . 
Changing to momentum space, Eq. dXT|) yields 

n ' = jf E E^ k > k 04 s ct k+qs ,c_ k , +qs ,c k , s , (A2) 

k,k',q s,s' 

where u(k,k') = v(k — k') due to translational symme- 
try of the crystal. Since we are interested in a situation 
with two sites per unit cell, we introduce two species of 
electron operators, 



„t 



-k+Qs 



a = 1, 
a = 2, 



(A3) 



where Q = (0,0, n) for a system as described in sec- 
tion IIIII and Q = (tt, tt) for the situation of sections IIVI 
andfVl respectively. Accordingly, we restrict the sum in 
Eq. (|A2[) to the two cases of q = and q = Q in the 
following. For the case q = we find 

n 'o = 4 EM k - kO[4 s 4_ ks , Cl _ k , s , Clk , s 



kk' 



C 2ks C 2-ks' C 2-k's' C 2k'. 
t J 



+ V{k - k' + Q)[cIk S C l-k S ' C 2-k' S ' C 2k' 



4k S C 2-k S ' C l-k' S ' C lk'J}- ( A4 ) 



For the other case q = Q, the interaction term can simi- 
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larly be written as 

^Q = 4EM k - k ')[ c L4- 



ks'^a-k's'nk', 



kk' 



C 2ks C l-ks' C l-k's' C 2k's] 



+ u(k-k' + Q)[c T 1 



^ks C 2-ks' C l-k's' C 2k'. 



+ C 2k S C l-k. 5 ' C 2-k' S ' C lk'J}- ( A5 ) 

At this point, we can distinguish the two cases of an 
interaction between sites belonging to the same sublattice 
and between sites on different sub-lattices. For the former 
case, i,j e A (B), we can use v(k + Q) = v(k) to write 
the above expressions as 



^o. Q = ^E^-( k - k >i 



ks L ^-ks'S-k's' L <5k's' 



k,k' 



with 



(A6) 



v% lS = v(k-k')[(r°U(r%}, (A7) 
«2s 7 * = ^k-kOKr 1 )^^ 1 )^]. (A8) 

Similarly, for the latter case, where i £ A(B) and j 6 
B(A), v(k + Q) = -v(k) yields 

v% l5 = v(k-k>){(T 3 U(T%}, (A9) 
«3w = v(k-k')[^) a ^T%]. (A10) 

In addition, the interaction can also be separated in a 
spin-singlet and a spin-triplet channel introducing Pauli 
matrices for the spin degrees of freedom, 



/ , C aks C £-ks' C 7-k's' C <5k's 



= 9 E A ^ 



1 W 



> C 7-k'S3 C <5k's 4 ' (AH) 



2 --°i s 2S3S4 u a ksi /3-ks 2 7-k's 3 L -(5k's4' 

S1...S4 

where 

A S1S2S3S4 = (c°) SlS2 (?°)U + (0.x* • «>U- (A12) 

Here, we have introduced ^° — ia y and <f = cri<7 y for 
simplicity of notation. 

The total interaction term now has the form 



'ks2 C 7-k's 3 C (5k's4- 



k,k' 



(A13) 

The interaction matrix element [V(k, k')]^^ 384 has an 
odd and an even part in k which depends on the result- 
ing sign of an interchange of the two first index pairs, 
(a/3, sis 2 ) «-» (/3a,s 2 si), 



WIV l r ' N ll s l s 2S3S4 _ ,,+ /\SlS 2 S3S4 



+ «kk'A!!^T/, (A14) 



where 

<k' = ^Kk-kO±^(k + k')). 
For the intra-sublattice interaction, these read 

ASXS2S3S4 _ ( 0\ I 0\f x 

xl^UirXs + ^UiT 1 )^} (A15) 

and 



AS1S2S3S4 



x[(t°U(t% + (t% p (t%], (A16) 
while for the inter-sublattice interaction, we find 

+ i(iT%Air% s m siS2 -&i Si (A17) 

and 



+ [(r 3 )Mr%ms lS2 -(<t)ls^ (A18) 

Unlike the case of a primitive unit cell, the momentum 
dependence is thus not only depending on the spin part 
of the interaction. 

As an example, we look in the following at the spe- 
cific example of stacked layers of Sec. IIIII The simplest 
non-trivial intra-sublattice interaction is between nearest 
neighbors, i.e., v(k— k') = V [cos(k x — k' x )+cos(k y — k' y )] = 
-<u(k-k' + Q). 

v ik' = {cos k x + cos ky) (cos k' x + cos k'y) 

+ ^-(cosk x — cosk' )(cosk x — cosfc| ') (A19) 



and 



v k k , = — U(sin k x sin k' x + sin k y sin k' y ). (A20) 

Note that for the cases of Sees. IIVI and |V1 the above 
functions correspond to the inter-sublattice interaction. 

For the nearest-neighbor inter-sublattice interaction, 
we find 



and 



w kk , = V cos k z cos k' z 



v kk , = — V sink z sin k' z . 



(A21) 



(A22) 
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Appendix B: Hopping matrix elements in systems 
like Fe-As-compounds 

The special structure of the FeAs layers in the iron- 
pnictidcs leads to a spin-orbit coupling with a differ- 
ent sign depending on the sublatticc. In this appendix, 
this spin-orbit coupling is derived for a simplified orbital 
structure, considering s-like orbitals for the Fe sites and 
p-type orbitals for the As ions, by focussing on only one 
sublattice [see Fig.UJa)] . To analyze the nearest- neighbor 
hopping - corresponding to a next-nearest-neighbor hop- 
ping in the full structure - it is easiest to rotate the crystal 
by 45 degrees and start with the As ions first lying on 
the bonds [see Fig.@|b)]. For this situation, the electrons 
can only hop from one Fe to the next in x (y) direction 
via a p x {py) orbital with hopping element t sp , 



fir, 



(x) 
sfi+x/2s 



t (») 
C isPi-x /2s 



c f v {v) 



cWi/j+h-c.: 



Hsl J i+y /2s ~ Hsfi-y/2s) "-^-J" ( B1 ) 

Assuming that the As-orbital's on-site energy differs 
from the energy of the Fe orbitals, Ea s = Ep c — A, we find 
for the nearest-neighbor-hopping integral in the effective 
one-band model 



' s\> 

A"' 



t' = (B2) 
The Hamiltonian in momentum space thus reads 

Hh op = ^ e hop c t sCks; (B3) 
k' 

where e{^? p = — 2i'(cos k' x + cosk' y ) with the new rotated 
axes k' x and k' y . Rotating the crystal back by 45 degrees 
to change to the old axes we find using k' x = (k x — k y ), 
k'y = (k x + k y ) and 

cos(k x ± k y ) = cos k x cos k y =p sin k x sin k y (B4) 

the usual (nnn) hopping energy e£ op = —At' cos k x cos k y . 

If the As ions are moved out of the plane, it becomes 
also possible to hop via a p z to a neighboring Fe site with 
hopping integral t sp . We therefore find the additional 
hoppings 



n = -t sP Y^sP% x/ 2s + 4sp¥ x/ 



2 s 



t («) 
C isP i 



i+y/2s ^isf l - y /2s> 1 " ,V "J- ( B5 ) 

We can now change to eigenfunctions of the As-site SOC 



with t = (t sp + it sp )/ y/2. 

Again reducing this to a single-band model by inte- 
grating out the As orbitals, we find in addition to the 
hopping Hamiltonian 



(B7) 



with 



<h3> s 



t (tsp ^sp) 



sp "spj ^2 _ \2 



(B8) 



a new SOC term, 



n' = J2 ( iac \s° V ss< c i+xs> - ia 4s a ss' 1 



iac\ s a x ss ,c l+ ~, - iacla^c.,, + h.c. ) (B9) 



with 



2t S pt S p\ 
A 2 - A 2 ■ 



In momentum space, this additional term reads 



k.s.s' 



(BIO) 



(Bll) 



where Ak = 2a(xsu\k v — y sin k x ). 
To transform this back, we use 

sm(k x ± k y ) = sin k x cos k y ± cos k x sin k y (B12) 

and also the rotated Pauli matrices, 



a* ^ ^(a*- a y), 



(B13) 
(B14) 



p) s , where the spin-quantization axis has to be orthog- Finally, we find the SOC Hamiltonian 
onal to the hopping direction to find 



^ = -E(* c U3 /2s + **44;i/2 



"Hsoc = E (^ k ' ^8*') 4 



k.s.s' 



-ks°ks'> 



(B15) 



t* C \sP\- x /2s * C lsPi-i/2s 
- (MsP^s ~ WAsV%\l2s 

(i*)*4p i % 28 + (i*)cLpU/J + h.c.) 



where now Ak = a(x sin fc^ cos k y — y cos k x sin k y ) . In 
a crystal with D^h symmetry, this term belongs to the 
irreducible representation B\ u . 



(B6) 
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